Optimal quasi-free approximation: 
reconstructing the spectrum from ground state energies 
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The sequence of ground state energy density at finite size, e l , provides mucii more information than usually 
believed. Having at disposal for short lattice sizes, we show how to re-construct an approximate quasi- 
particle dispersion for any interacting model. The accuracy of this method relies on the best possible quasi-free 
approximation to the model, consistent with the observed values of the energy bl. We also provide a simple 
criterion to assess whether such a quasi-free approximation is valid. As a side effect, our method is able to 
assess whether the nature of the quasi-particles is fermionic or bosonic together with the effective boundary 
conditions of the model. When applied to the spin- 1/2 Heisenberg model, the method produces a band of Fermi 
quasi-particles very close to the exact one of des Cloizeaux and Pearson. The method is further tested on a spin- 
1/2 Heisenberg model with explicit dimerization and on a spin-1 chain with single ion anisotropy. A connection 
with the Riemann Hypothesis is also pointed out. 



I. INTRODUCTION 

The way finite-size, thermodynamic, averages approach 
their infinite size limit, encodes a wealth of precious infor- 
mation. Take for example a quantum system at zero tem- 
perature, and consider the ground state energy density at fi- 
nite size bl, L being the linear size of the system. A piece 
of common knowledge in condensed matter theory -not yet 
a theorem though- asserts that, for a massive (i.e. gapped) 
system with translation and Left-Right invariance, the ap- 
proach to the thermodynamic limit is exponentially fast, that 
is — ?> Coo + O (e^'^/^^) where the decay rate is some- 
how related to the correlation length of the system £^c- On 
the contrary, in a critical theory, the correlation length is for- 
mally infinity and the approach is expected to be of algebraic 
type. Moreover for conformal field theories (CFT) in (1 + 1) 
dimensions and periodic boundary conditions, it is known that 
gl = Goo — (7r/6) cv/L'^ + o (L~^)i2 where v is the speed 
of elementary excitations and c is the all-important central 
charge of the CFT. From these argument it seems clear that 
it should be possible, by simply looking at the sequence e^, 
to discern whether the theory is critical or gapped. In fact, 
even though the standard method to locate critical point is that 
of looking at the gap closure, recent studies have shown that 
accurate methods to locate critical points are available, which 
resort solely to the computation of ground state averages such 
as the total energy and the perturbation^. 

In this article we will show that the sequence e l provides in 
fact much more information. We will show that, by knowing 
cl for few short sizes L, it is possible to reconstruct approxi- 
mately the one particle dispersion of the theory. Perhaps most 
importantly, analyzing the sequence e^, one is able to assess 
whether the quasi-particles of the theory have fermionic or 
bosonic character together with the effective boundary condi- 
tions of the model. 

If the Hamiltonian H under consideration admits a quasi- 
free representation (i.e. the Hamiltonian can be expressed as 
a quadratic form in Fermi or Bose operators), our algorithm 
produces a one-particle dispersion rapidly converging to the 
exact one, as we increase the number of available ground state 
energies. In practice having 10 energy data produces a result 



indistinguishable from the exact one. 

Conversely, if H is truly interacting, we obtain an approx- 
imate one-particle dispersion corresponding to a quasi-free 
model H. The Hamiltonian H so obtained, is the simplest 
quasi-free model whose ground state energies are the same as 
those observed for H. 

For quasi-free models the mathematical problem is related 
to that of reconstructing a function given some of its partial 
Riemann sums. As we will show, assessing the convergence 
speed of this method is a problem connected to the Riemann 
hypothesis. Let us then begin by considering the problem of 
convergence of Riemann sums. 



II. CONVERGENCE OF RIEMANN SUMS 

As we will show in greater detail in section IIVI the finite- 
size, ground state energy of any, translation invariant, quasi- 
free system, consisting either of bosons or of fermions, is re- 
lated to a Riemann sum of the following form 
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Here L is precisely the size of the one-dimensional latticai^ 
and / a suitable function related to the one-particle dispersion. 
The points A;„ = (27rn + ■&) /L are quasimomenta which de- 
fine the Brillouin zone BZ = [0, 27r)''' for finite L, and the 
angle d G [0, tt] defines general, twisted boundary condi- 
tions (TBC). In terms of canonical operator (either bosons or 
fermions) Cx, TBC means Cx+l = e'^'ca;. The angle d allows 
to interpolate continuously between periodic (PBC) i9 = 
and antiperiodic (ABC) d = n boundary conditions. 

Since the function / is defined on [0, 27r) it is useful to write 
it as a Fourier series: 

oo „2ir 

/(fc)= c„e"^ c„ = — / e-^«"/(Ode (2) 

In order to obtain an explicit formula for Eq. ([U we would 
like to apply the operator Sl to the exponentials appearing in 
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Eq. ^ (i.e. swap Sl with the sum). To this end we need to 
assume some regularity in the convergence of the sum (|2| such 
as absolute convergence. Necessary conditions for absolute 
convergence of Eq. (|2|l are known; for instance it is enough to 
have /, a-Holder continuous with a > 1/2, or / of bounded 
variation and a > 0. On physical grounds we can assume 
such conditions will be satisfied, they are indeed satisfied in all 
examples encountered. Now, assuming absolute convergence 
of the series in Eq. we can apply the operator Sl to the 
exponentials obtaining 
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This equation is zero for all natural n not multiple of L, 
whereas for n — IL with I e Z, equals e*'"'. Therefore we 



can compactly write Sl (e™'^) = X^z^--- 
bining this result with Eq. (|2| we obtain 
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SL{f)-co = ^e'^L+e-^'^'c-zL 
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where the last line follows from the reality of /. 

Assume now that the function / is 27r-periodic and ana- 
lytic. Then there exist a strip in the complex plane A = 
[0, 2tt] X cr, a] such that / can be extended to an analytic, 
bounded function on A. In this case Fourier coefficients can 
be computed in the complex plane z, integrating on a hori- 
zontal line shifted by an imaginary amount — icr. The vertical 
contributions cancel each other because of periodicity and we 
can write 
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7(z)dz. 



(6) 



At this point the geometric series of the exponential con- 
verges absolutely and we obtain 



/ closest to the real axis, let it be zq. The imaginary part 
= Im (zq) = max a gives the decay rate and from Eq. ([8]) 
we obtain \Sl (/) - co| < O (e--^/«^). 

As we will discuss in greater detail in section [TVl for quasi- 
free systems (that is, Hamiltonians quadratic in Fermi or Bose 
canonical operators), the ground state energy is exactly given 
by an expression of the form of Eq. ([T]i where the function / is 
proportional to the one-particle dispersion. The result above 
shows that if the one particle band is a positive, analytic, pe- 
riodic function, the ground state energy density approaches 
its limit exponentially fast. The positivity requirement on the 
band implies that the system has a gap. So we reach the con- 
clusion that for gapped theories with analytic dispersion, the 
energy density decays at least exponentailly with the size -at 
least for the case of quasi-free systems. Since, from scaling 
hypothesis, the correlation length is the only length-scale of 
the system, we expect the decay rate to be related to the 
correlation length. 

Equation (|7]) cannot be used when the function is not an- 
alytic on [0, 2tt], in this case the difference 5*^ (/) — cq can 
be estimated using Euler-Maclaurin formula. In any case, 
for critical theories where the dispersion vanishes linearly as 
Lj (fc) ~ v\k — kpl'^ (here z is the dynamical exponent and 
kp the Fermi momentum) scaling theory predicts an algebraic 
approach of the form 5*^ (cj) — cq = O (L^^^^) (see e.g.'*). 
An example in this class is given by the function |cos(/c)|, 
which is continuous but not differentiable. An explicit calcu- 
lation shows in fact Sl (|cos {k)\) — 2/tt ~ O {L^^) consis- 
tent with z = 1 . To complete the scenario we should mention 
another possibility which gives rise to an algebraic approach 
to the thermodynamic limit. Namely the function itself might 
have a jump at some point as it happens, for instance, to the 
function k. In fact and 2tt should be identified in the Bril- 
louin zone and this function has a jump of value 2tt. In fact in 
this case we have Sl (k) = tt — ir/L. This situation typically 
takes place in physical systems when Left-Right symmetry is 
explicitly broken. 

Our aim here is to reconstruct / given Sl (/), that is, we 
would like somehow to invert Eq. To this end it is better 
to do a step back. 



Sl if) - Co = Re 



/ (z) dz 
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III. SOLUTION OF THE INVERSE PROBLEM 



This formula is remarkable; it allows to show that, for periodic 
analytic functions, the convergence of partial Riemann sum to 
the limiting integral, is at least exponentially fast in the size 
L. More precisely, since 1/ |e'(^-^-'') - l| < 1/ (e'"^ - l) 
for all z in A and — Q,tt, using dTjl we obtain 



|5l(/)-co| <2max|/(z)| 
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This shows that, for analytic periodic functions /, the par- 
tial sum Sl (/) converges at least exponentially fast in L 
to its limiting value. Moreover Eq. dHJ gives a simple way 
to compute the decay rate. All we have to do is take cr as 
large as possible, i.e. find the largest strip A where the func- 
tion is analytic. This amounts to look for the singularity of 



Even in the most favorable case, we can hope to have at 
disposal only a finite number of partial sum S,n if), say the 
first L (obtained, for instance by Lanczos diagonalization, see 
sec. IIVI ). This means we can only reconstruct / through L 
coefficients. A sensible choice, and in some sense the best 
one, is that of reconstructing the first L Fourier coefficients of 
/■ 

In the rest of the article we will concentrate on the phys- 
ically most important periodic and anti-periodic EC's for 
which i9 = 0, TT. In these cases the partial Riemann sums 
are identically zero for / odd, i.e. S^j^^ (/) — for / (fc) — 
— / (— fc) [~ — / {2tt — k)]. Clearly, having at disposal the 
numbers sf ^ (/) we can only hope to re-construct the even 
part of /. Fortunately, on physical grounds we can suppose 
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that the function / will indeed be even in k. This corresponds 
to Left-Right symmetry on top of translational symmetry and 
it is a reasonable assumptions valid in most physical situa- 
tions. A breaking of Left-Right symmetry happens, for in- 
stance, when inserting a periodic spin chain in a magnetic 
field of flux dhc/e^. Such a problem can usually be reformu- 
lated into the same problem without magnetic field and TBC 
specified by the angle -d. All in all considering i? = 0, vr is 
consistent with assuming / even. 

Assume then that the function / is even so that its Fourier 
series has only cosine coefficients. For PBC and ABC, e"' 
is real and, defining Rl (/) — Sl (/) — cq, we can re-write 
Eq. (HI in matrix notation as il = Qa, where boldface indi- 
cates column vector and the matrix Q has components given 
by 



The above sum extends over all positive divisors of m 
and is usually denoted by X]j|m (r^^d j divides m) in 
the mathematical literature. Equation ( fTsT l provides a re- 
cursive solution for h (m): b (1) = e^''' while h (j) = 

-e'^'E^ ^li e*''J'/™&(TO) for j > 2. Since Eq. O is 

j/meN 

independent of L, it implies that [^'^^] is the first L x L 
sub-matrix on the diagonal of [t/'*^'] ^ for any L < M. In 

particular any [^'^-'j ^ can be obtained from the infinite case 
Q^^ and its entries do not depend on L. Now the inverse for- 
mula a'^' [^'^'] i?^^', in components reads 



[L/k] 

4"^' (/)=EM")i?nfe(/): 



(14) 



1=1 



m,lM ■ 



(9) 



The components of the vector a are the cosine Fourier coeffi- 
cients 



2Rec„ 



1 



cos(en)/(e)de 



(10) 



Note that both Q and a are real. To solve for the first L Fourier 
coefficients we truncate the equation to the first L terms to 
obtain = g('^)a('^), where J?^^) = (i?i, i?2, . . . , RlV , 

is the L X L matrix wil 
< L. For example, for 



a^^") — (ai, a2, . . . , a^)^ and Q^'"'' is the L x L matrix with 
entries given by Eq. (|9]l for AI, m 
L = 6 , the matrix Q''^^ reads 
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(11) 



with q 



±1. Since det (g^^)) = q^ =^ the ma- 
trix is invertible and defines a bijection. It is then possible 
to obtain the first L Fourier components approximately via 

a(^) = [S^^-*] i?*-^-*. One can show that the inverse matrix 
has the same structure as Q^^^ in the sense that [^/^'^^j ^ has 
non zero entries in the same places as Q'^-'^K More precisely 
b {j /i) if j/i G N and zero otherwise, for 



some numeric function b (n). Using a similar notation as be- 
fore, this means 



T(i) 



[L/i] 
'■j (=1 



(12) 



Imposing [C/'^^-'] ^ 5'^-' = I we get the equation defining 

b (m): 



E 

jn = l 



(13) 



where we explicitly indicated the dependence on the function 
/. The superscript (L) indicates that the Fourier coefficients 
are obtained only approximately, but with increasing precision 
the larger the L. The result of the reconstruction is optimal in 
the sense that it produces the unique trigonometric polynomial 
of degree L consistent with the "observed" data Sl- In the 
limit L ^ oo - that means we know Rl (/) for arbitrary L - 

we can reproduce the function exactly and a^^^ (/) ak (/). 
To be mathematically precise, this last assertion is satisfied 
provided / is not too pathological. On physical ground we 
can safely discard such pathological cases. On the contrary 
Eq. ( fT4] i can be used to obtain information on the function 
/. For what discussed in the previous section, the only po- 
tentially dangerous case is that of a critical point. Since 
physically we will identify Rl (/) with the finite size en- 
ergy density, scaling arguments predict that, at a critical point, 
Rl if) ~ L~^~^ where z is the dynamical critical exponent. 
Now one can use Eq. (fT4l l to obtain k^^^'- which gives 
rise to a function with absolutely convergent Fourier series. 

The PBC case ?9 = provides some interesting connec- 
tions to number theory. In this case in fact Eq. (fTsT l be- 
comes the equation defining the Mobius function /Lt (m)^, 
i.e. Eq. (fT3] l is solved by b (n) — ji [n). Note also that, for 
q = 1, Q'^^^ is very similar to the Redheffer matrix TZ'^^^ 
known in number theory^. The matrix TZ^^^ is the same as 
(y(^) except for the first column which is made of one. The 
importance of the Redheffer matrix originates from the fact 
the det7^(^) = M(L) := where M (L) is the 

Mertens function. The statement M (L) = O (L^/2+^) is 
equivalent to the Riemann hypothesis. 



A. Convergence rate 

Since we are assuming the function we seek is even, Eq. (|2|l 
becomes / (fc) = cq + cos [nk). If we have access 

to the partial Riemann sums up to L (and the limiting value 
Co), we can reconstruct an approximate function given by 



/l (fc) = CO + E "i"^^ cos(nfc). 



(15) 
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with ai^^ given by equation ( IT4] |. We can now ask how fast 
this method allows to reproduce the function /. The question 
of convergence of Fourier series has engaged mathematicians 
for centuries. The problem of convergence of the series ( fTSl l 
is likely to be more complex. Here we will content to give 
some arguments for the physically important cases related to 
massive and critical theory. 

To be specific we will consider PBC. We can safely assume 
/ to be square summable. The reconstructed function /l is 
obviously also square summable, being a trigonometric poly- 
nomial. The ([0, 2tt]) square distance reads 



2v 

\.fL-.f\\l = 7rY.{c 



i=L+l 



(16) 



Consider first the case where / is periodic and analytic in 
[0, 27r]. In this case a saddle point argument shows that the 
Fourier coefficients decay exponentially c„ — 0(6^"/^^^). 
Moreover the correlation length S^p is precisely the same as 
that appearing in Sec. [Ill = S,e- Hence the second sum in 
Eq. ( fT6b is of the order of e~^^/^^. To estimate the first sum 
consider 



oo 



m=[L/n] + l 



b (m) R„ 



For what discussed in Sec. |II1 = O 

PBC \b{m)\ = |M("^)| < 1 so that 
^m=lL/n]+i \Rnm\- This implies that, for sufficiently large 



(17) 



For 



< 



l^e ■^/''^ ) , so that, in turn, the first sum in 

Eq. ( fT6b is bounded by Le^^^/^^. All in all, if / is periodic 
and analytic in [0, 27r], is exponentially close to /, in sense 
that||/i-/||^<0(Le-2V«-). 

Consider now the critical case. From a physical point of 
view, a critical theory with dynamical exponent z corresponds 
to excitations vanishing as / (k) ^ k^, with z > 0. A scaling 
argument now implies that the Fourier coefficients of such a 
function scale as c„ ~ n^^^ 
for large enough n |c„| < An 
Let us first consider the second sum in Eq. (fT6b 



More precisely assume that, 
with A positive constant. 



n=L+l 



< A- 



n=L+l 



We can use the asymptotic behavior 



(a - 1) ' /M"^ 



to estimate 



(18) 



< 



A^{l + 2z) ^ /L^+^''. To obtain the behavior of the 
first sum in Eq. ( fTSl l we first look at Eq. Q and obtain, for L 
large enough \Rl \ < 2AL-'^-% (1 + z). Then 



< 



2AC{l + z) 



,1 + 2 



E 

= [L/n\ 



Using the same estimate as before we obtain 



2AC(l + z)z-V(nL^) and so ^ 



L 



(19) 

< 
< 



A'/{L^^) < ^'7(L22). So it seems that the con- 

vergence rate of — fW^ is of the order of L^^^. Here we 
would like to conjecture that the first sum in Eq. ( fTSl l actually 
introduces corrections that are roughly of the same order as 
the second, i.e. — f\\l ^ L~^~^^. The argument is based 
on the Riemann hypothesis. Consider again — a„ in the 
PBC case. For L sufficiently large we have 



(to) to ^ ^ . 

m=[L/n] + l 



(20) 



An equivalent statement to the Riemann hypothesis is that the 
Mertens function AI (x) := Em=i A* ('^) satisfies M (x) = 
O for any positive e. Using partial summation^' we 

can estimate Em=Af 1^ ("^) = O (A/^/^^"+'^) for a > 
1/2. Plugging this estimate in Eq. ( f20l i we obtain 



„l/2+e 



(21) 



-l-2z+e 



from which roughly E^=i ~ O")^ ^ L 

To summarize, in the case of analytic functions, relevant to 
massive Left-Right symmetric models, we expect exponential 
convergence speed with rate given by the correlation length. 
For critical theories with dynamical exponent z, we expect 
an algebraic convergence speed. Resorting to the Riemann 
Hypothesis we conjecture the convergence speed to be of the 
order of L-l-2^+^ 

To visualize better how fast the method works, consider an 
example taken from physics where the function to reconstruct 
is given by 



/(fc)- J^sin(fc/2)' 



(22) 



Such a function is a periodic generalization of a relativistic 
dispersion and it describes, exactly or approximately, the one- 
particle dispersion of many one dimensional systems. When 
TO 7^ the function is analytic in [0, 27r] and the convergence 
is exponentially fast. The case to = can serve to model 
critical theories with dynamical exponent z = 1. In this case 
the convergence in the L^-norm, is algebraic and we just con- 
jectured that the rate is of the order of L^^. This behavior is 
confirmed in figure [T] which shows that, for the massive case, 

ln(||/L 

TO = we have roughly In — f\\2j 

The result of the reconstruction for the function in Eq. 
is instead shown in figure|2|for different masses to, using PBC 
and as little as the first ten Riemann sums. Results for ABC 
are very similar. It is notable the very good agreement even in 
the massless case. 



IV. METHOD AND APPLICATIONS 

The methods discussed so far are readily applicable to 
translation invariant quasi-free systems consisting either of 



/II2) is approximately linear with L, while for 

< — 31nL + const., 
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Figure 1: Convergence speed in the ([0, 2tv]) norm for the func- 
tion / (k) — \J sin (fc/2)^ + m^. The reconstructed function 
/l (fc) is given by Eqns. l |15t and ( I14t and the ([0, 27r]) distance 
1/ — /lIIj is then evaluated for various L. Top (bottom) picture 
refers to massive m 7^ (massless ra = Q) case. In the massive 
case the approach expected is of exponential type. The dashed line is 
given by const, x L~''^ e~'^'"^^^ , with ^^j^ = 2arcsinh (m) which 
can be obtained estimating the Fourier coefficients of / (k) with a 
saddle point argument. In the masless case the approach is expected 
to be of the form l/L'^^'^ (see main text). The dashed line in the 
lower panel gives the behavior. 



and the mass in the top panel is m ; 



The method is tested for PBC, 
0.1. 



fermions or bosons. By quasi-free systems we mean here 
Hamiltonians that can be expressed as quadratic forms in Bose 
or Fermi operators. In such cases in fact the ground state en- 
ergy is precisely given by a partial Riemann sum. For ex- 
ample, in the notation of^ a quasi-free fermionic model has 



the form H = c]A 



1/2 



Di- 



agonalization brings it to = Yl,k '^kvlvk + T, where the 
band utk can be chosen positive and the constant T is given 
by 2r trA — J^k^k- For translation invariant systems, 
with PBC or ABC, the label fc is a (quasi-) momentum quan- 
tized according to fc = (2?™ + d) /L, n = 0,1, . . . , L ~ 1. 
Defining the "filling fraction" ly = 1 — tr {A) / ^k the 
Hamiltonian takes the form 



(23) 



If the model consists of n species of non-interacting colors, 
i.e. Ci Cj Q, a = l,...,n, simply replace ly nv. 
All the methods presented so far can now be applied con- 
sidering that the ground state energy density of Hamiltonian 
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Figure 2: Results of the reconstruction for a prototypical function 

/ (fc) = J \J sin (fc/2)^ + irfi, using as little as ten Riemann sums. 
Upper panel refers to m = 0.1 while lower panel to m = 0. Con- 
tinuous line is the exact function, while the dashed line is the re- 
constructed function (fc) using PBC. Since the functions are even 
only half of the Brillouin zone is shown. Note that for m = 0.1, 
/l (fc) is indistinguishable from the exact function at this scale. Re- 
sults for ABC produce very similar plots. 



is precisely given by cl = El/L = - {v/2)Y,k^k- 
Now the point is that in many physically interesting situa- 
tions the "filling fraction" is known in advance. In fact in 
absence of (magnetic or electric) fields generally v = \ (or 
V — n for n non-interacting species) since trA = 0. This 
means that the v in Eq. ( l23T l is independent of L, and the 
ground state energy density is precisely given by a Riemann 
sum: = — {v/2)Sl {(^k)- Similar considerations hold for 
bosonic quadratic theory with the important difference that 
now cl = + [y 12) Sl {(^k), due to the commutation relation. 

We can now argue that any interacting model admits some 
sort of quasi-free approximation. At this level of approxima- 
tion, the Hamiltonian is quadratic and we can apply all the 
reasoning presented above. Our method gives a way to ob- 
tain a one particle dispersion knowing the ground state energy 
for some lattice sizes. The dispersion obtained is optimal in 
the sense that it is the unique trigonometric polynomial of de- 
gree M (where AI is the number of energy data) consistent 
with the observed values of the energy. This method bears 
some similarity with the Hartree-Fock method largely used 
for ab-initio calculation of molecular systems. The Hartree- 
Fock method, for a given size L, gives the optimal quasi-free 
state that minimizes the energy. The method proposed here 
instead, taking AI ground state energies as input, gives an op- 
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timal quasi-free system (identified with its one-particle disper- 
sion), in the sense that its ground state energies are precisely 
the M observed value. 

To specify completely the problem one has to assume the 
character of the quasi-free approximation, i.e. whether the 
model consists of Bosons or Fermions together with the ef- 
fective boundary conditions. In practice we have to chose 
if the ground state energy densities are given by — 
e (zy/2) 5*1, (wfc) with e — ±.1 and moments specified by 
1? = 0, TT (PBC or ABC). This choice can be straightfor- 
ward if the model under consideration consists of Bosons or 
Fermions, but in case of spin models, the character of the 
effective, quasi-free model is less clear. According to the 
choices e = ±1 and d = 0,7r we have therefore 4 possi- 
bilities. However the requirement that the reconstructed band 
must be positive fixes in practice only two combinations. This 
is an important result on its own: simply looking at the se- 
quence of ground state energies, one is able to assess whether 
quasi-particles are Fermions or Bosons with ABC or PBC. 

As for any approximation method, it would be derisable to 
have a simple criterion to assess whether a quasi-free approx- 
imation is feasible in the first place. Such a criterion can be 
given. In fact for exactly quasi free models, the ground state 
state energies satisfy 



where the superscript (0) , (tt) refers to PBC, ABC respec- 
tively. So, having finite-size energies for PBC and ABC, we 
can simply verify the possibility of an effective quasi-free de- 
scription by checking how well Eq. ( l24l i is satisfied. 

For what we have said in section IIII Al the procedure of 
reconstructing a function given its M partial Riemann sums, 
rapidly converges upon increasing M even in the the critical 
case (the worst scenario), so that one can effectively limit one- 
self to short lattice sizes. 

For the reader's sake, let us sketch here the relevant steps of 
the algorithm: 

• Obtain a set of M ground state energies of the system 
by exact diagonalization of short lattices, say sizes up 
to L = 10 - 20. If both PBC and ABC energies are 
available one can check the feasibility of a quasi-free 
approximation by checking how well Eq. ( |24] | is satis- 
fied. 

• Assume effective PBC/ABC and Bosons/Fermions 
which correspond to assume for the ground state energy 
density — e {v/2) Sl (wfc) with e = ±1 and mo- 
ments specified by — 0,7r. An approximate disper- 
sion is then given by Eq. ( fTSl l with f (k) — e {v/2) uJk 
and coefficients specified by Eq. ( fT4l i. The requirement 
u! {k) > will fix two cases out of the four possibilities. 

• One should also fix the filling fraction v. If this can be 
simple when the model is originally given in terms of 
Bosons or Fermions, in general one must be guided by 
physical intuition. Referring to the example that will 
be discussed in the following sections, it is natural to 
expect V — 1 for the pure spin- 1/2 Heisenberg model. 




Figure 3: (Color online) Results of our procedures for the spin-1/2 
Heisenberg model. The continuous upper curve (blue) is the ex- 
act dispersion uj (k) —{n/2) J |sin (fc)|, while the continuous lower 
curve (green) is the result of spin-wave approximation ujsw (fc) — 
J |sm (fc)|. Using only finite size energies up to L = 24 we obtained 
upper and lower dashed curve. The lower curve is obtained assuming 
a bosonic dispersion with effective PBC and can be identified with 
the spin-wave approximation. The upper dashed curve assumes a dis- 
persion with fermionic character and effective ABC. The very good 
agreement with the exact result tells us that there exist a quadratic 
Hamiltonian approximating Heisenberg one very precisely. 

a triplet of excitations for its dimerized version (ly = 3) 
and a doublet {ly ~ 2) of excitations for the spin-1 X — D 
model in the large-Z? phase. 

Let us now illustrate how the method works on the hand of a 
few concrete, yet prototypical examples. 

A. Spin 1/2 Heisenberg model 

Take the Heisenberg antiferromagnetic (J > 0) chain: 

L 

i7= J^S,;-S,;+i. (25) 

i=l 

Si are spin-1/2 operators at site i and PBC (S^+i = Si) are 
used. From the exact solution we know that the infinite size 
ground state energy is Boo/"/ = 1/4— In 2**, whereas the quasi- 
particle dispersion is given by uj (k) = (7r/2) J |sin (fc)|^. 
When we have to evaluate the energy at finite size we immedi- 
ately face a problem. When L is even the ground state belongs 
to the total spin S ~ sector and is unique'". This is the kind 
of ground state we "expect" from this model. On the contrary, 
for L odd there are two degenerate spin-1/2 ground states. The 
ground state energies for L odd have a completely different 
character and our intuition suggests us to discard them. As a 
consequence we have access only to for even L. However, 
from the exact solution we know that the dispersion is peri- 
odic with period halved i.e. tt. Hence it has only even Fourier 
(cosine) coefficients. From Eq. ( fT4b we see that with even- 
size energies we can re-construct even Fourier coefficients. In 
this case the two facts are consistent: e l only for L even ^ 
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even Fourier coefficients. The procedure is as follows. First, 
we diagonalize exactly Hamiltonian ( |25] l with say a Lanczos 
algorithm. In few minutes of a small laptop computer, we 
obtained ground state energies for lattices of even size up to 
i = 24 . Separately we estimate the infinite size ground state 
energy which in this case is Coo = «/(l/4 — In 2). Having 
collected the numbers for L even, we can use Eq. (fl4] i 
to obtain a one-particle dispersion. To specify the problem 
completely we have to make few further assumptions. First 
we have to fix boundary conditions. Even if we have PBC 
for the spins different EC's can be induced in the effective 
quasi-free model. Indeed using the Jordan- Wigner transfor- 
mation, model dZSl l can be exactly mapped to a model of inter- 
acting spinless fermion with parity dependent boundary con- 
ditions (see for example" for a discussion on these emerg- 
ing BCs). Since the ground state is a singlet it belongs to 
the parity one sector, where BCs for the fermions are anti- 
periodic. So, to be more general, we consider equation (fl4] i 
for (7 = 1,-1 which corresponds to effective PBC or ABC 
respectively. On physical grounds^- we fix the filling fraction 
to ly — 1. This is enough to obtain a q dependent function 
f (k) — (e/2)cjfc. To specify completely the dispersion we 
must still decide whether the effective quasi-particles are ei- 
ther fermions (e = — 1 ) or bosons (e = +1). The four possi- 
ble cases corresponding to q = ±1 and e = ±1 are reduced 
to two by imposing positivity of the band. The result is that 
assuming effective PBC quasi-particles are Bosons, while as- 
suming ABC quasi-particles must be Fermions. 



The results of the procedure, using only PBC ground state 
energies up to L — 24, are shown in Fig. (O. The bosonic 
dispersion with effective PBC can be identified with the spin- 
wave dispersion ujsw (k) — J|sin(A:)|, obtained with the 
spin-wave approximation of Anderson and Kubc'^ '^ [see 
lower curves in Fig. (O]. Instead, the Fermionic disper- 
sion with effective ABC is very close to the exact one of 
des Cloizeaux and Pearson. The excellent agreement of this 
dispersion with the exact one, indicates that a good descrip- 
tion (as long as short range quantities are concerned) of the 
Heisenberg model can be given in terms of an effective quasi- 
free fermionic Hamiltonian with ABC. 



Although the exact one-particle dispersion of the spin- 1/2 
Heisenberg model could be reproduced with high precision, 
this example also shows a limitation of our method. Namely 
to use our method we need either ground state energies for 
general L, or if we only have access to even sizes L we can 
only reconstruct even Fourier coefficients. These limitations 
disappears if we consider dimerized models where we expect 
the dispersion to be 7r-periodic. This is because an even func- 
tion of period tt has only even (cosine) Fourier coefficients. 
Knowing finite size energies for even sizes L is enough to 
reconstruct -within a certain approximation- the whole one- 
particle dispersion. 




Figure 4: One particle dispersion of model l i26b ai 5 = 0.048. The 
full dots are exact diagonalization data of^^ obtained for a chain of 
L = 28 sites (reproduction with permission of the authors). Our 
method is tested using only finite size energies up to L — 12. 
The dashed curve is the result of our method assuming a triplet of 
fermions with effective ABC, while the continuous curve assumes a 
triplet of bosons with effective PBC. 



B. Dimerized spin- 1/2 chain 

Consider then a spin- 1/2 Heisenberg model with an explicit 
dimerization of the exchange coupling: 



H 



1=1 



s? • s 



i+1 



(26) 



This model has been extensively used to characterize a vari- 
ety of spin-Peierls compounds. The presence of the dimeriza- 
tion has the effect of halving the Brillouin zone and so, cor- 
respondingly, the one-particle dispersion should have period 
TT. That this is indeed the case is confirmed by many numer- 
ical simulation''*. Then we can safely use even size energies 
to reconstruct the even Fourier coefficients of the dispersion. 
Moreover, a non-zero 6 has also the effect of opening a mass 
gap. As discussed in section IIII AI the convergence rate of 
our method is expected to be extremely fast in this case. With 
the aim of showing the usefulness of the method, we consider 
very short length. Using only finite size energies at even sizes 
from i = 2 to L = 12 (i.e. only six numbers!) we obtain 
the dispersion shown in Fig.|4] The results are then compared 
with those obtained via much more powerful diagonalization 
of refJi performed on a chain of i = 28 sites. 



C. Spin-1 model with single ion anisotropy 

As we have shown, our methods can be successfully applied 
to spin- 1/2 chains only when the dispersion is even and of pe- 
riod TT. This is the case for the pure Heisenberg model and for 
dimerized models as the one of Eq. ( l26b . What about spin- 
1 chains? For PBC and even size the theorem by Lieb and 
Mattis'" tells us that the ground state of a generic antiferro- 
magnetic Heisenberg model belongs to the total spin zero sec- 
tor and is unique. For odd sizes an antiferromagnet with PBC 
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Figure 5: One particle dispersion for model Eq. ( I27t in the large-D 
phase (D — 7.4). Only the part of the function corresponding to 
Fourier cosine coefficients c„ n > 1 are shown. The continuous 
curve is the result of third order perturbation theory Eq. ( I28t without 
the cos (k) term. Perturbation theory is basically exact for such value 
of D according tc'"". Our procedure is tested using finite size energies 
(obtained by Lanczos method) of short chains from L = 2, 3, . . .up 
to L = 12 and later discarding the cos (fc) term from the disper- 
sion. The dashed curve is our result assuming the dispersion is a 
doublet of bosonic particle and effective PBC. Assuming a doublet 
of Fermi particles and ABC one obtains similar accuracy although 
the L'^{[0,Tv]) norm distance with dispersion l l28t favors the bosons 
(distance 0.049 against 0.078). 



is frustrated and the theorem does not apply. However we have 
numerically verified that also for odd sizes the ground state 
belongs to the total spin zero sector (this is consistent with 
the VBS description and with the fact that every spin-1 can be 
thought of a symmetric combination of two spin- 1/2, and so 
any chain contains an even number of spin- 1/2). This suggests 
that we could use ground state energies both for even and odd 
sizes and re-construct completely the dispersion. However 
there is still a problem with this approach. The ground state 
energy for the single site problem is not clearly defined. If the 
model admits a quasi-free approximation, the L = 1 ground 
state energy is given hy Ei = ± (y jl) S\ (uk) = ± (i^/2) 
(plus or minus refers to Bosons or Fermion respectively). Us- 
ing the inversion Eq. ( fT4] i we see that, 5*1 (/) enters only in 
the definition of the first Fourier coefficient ai. So a missing 
Si (/) allows to specify the function up to an additive cos (fc) 
term. This term could be fixed by other means, such as ob- 
taining the value of the dispersion at a given momentum. Let 
us analyze a concrete example. 

Consider the spin- 1 model with single ion anisotropy 



H = J 



L 

E 

=1 



s, • s 



4+1 



(27) 



where Si are now spin-1 operator. In order to test our method 
we consider the model for large D where perturbation the- 
ory is applicable and an analytic expression for the disper- 
sion is available. When D = oo the ground state is given 
by |0,0, ...,0), and excitations form a doublet of degen- 
erate states with the spin at one site flipped to +1 or —1, 
i.e. |0, 0, . . . , ±1, . . . , 0). A finite large D, removes transla- 



tion degeneracy and one obtains a doubly degenerate band 
uj{k). This picture remains valid in the whole, so-called, 
large-D phase, which is separated by the Haldane phase 
roughly at Dc — 1 (see'^ for details on the phase diagram). 
In the large-Z? phase one can use perturbation theory to obtain 
the doubly degenerate dispersion. A third order calculation 
has been performed'^ with the result 



^=D + 2cosik) + ^ 



1 + 2sm{ky 



1 



2 sin (kf - -(1 + 8sin(fc)^)cos(fc) 



O D- 



(28) 



We re-write the dispersion as 



uj{k) 
J 



D2 



2 

cos {2k) 



1 



2 1 



cos (fc) 



£12 



COS (3fc) + O {D- 



(29) 



in order to make clear the Fourier (cosine) coefficients of 
the dispersion. Using a Lanczos algorithm we computed the 
ground state energy of the model (|27l) for L 2, 3, . . . , 12. 
For what we have said, using the inversion Eq. ( fT4b we can 
reconstruct the band up to a cosine term. In figure|5]we show 
the result for the reconstructed band compared to the disper- 
sion Eq. (|29] | both without the cos (fc) term and the agreement 
is excellent. As noticed previously, the cos (fc) term can be 
fixed by other means. 



V. CONCLUSIONS 

In this article we showed that a lot more information than 
currently believed, is encoded in the ground state energy den- 
sity at finite size e^. In particular we provided a method able 
to reconstruct an approximate one-particle dispersion for any 
one-dimensional quantum system, given some finite size nu- 
merical data {e^ }. The dispersion reconstructed with this pro- 
cedure is optimal in the sense that it is the unique trigonomet- 
ric polynomial of degree L,„ax (^max being the number of 
energy data) consistent with the observed data {e^}. Equiv- 
alently the method produces a quasi-free Hamiltonian which 
has the same ground state energy densities as the observed 
values {ei,}. This method is exact if the model has some 
sort of quasi-free representation, and it converges very rapidly 
increasing imax so that very few data are sufficient (using 
10 energy data gives already very good results). We also 
provided a simple criterion to assess whether such a quasi- 
free approximation is feasible in the first place. As a side 
effect, simply looking at the sequence {cl} this method is 
able to assess whether effective quasiparticles are either boson 
or fermions with effective periodic or anti-periodic boundary 
conditions. 

Since the Casimir force is specified (up to a constant) by 
the energies {bl}, from a physical point of view the proce- 
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dure presented consists on reconstructing the one-particle dis- 
persion given the Casimir force. 

Further developments in this direction include the possibil- 
ity of extending these ideas to higher dimension and testing 
the procedure on other strongly correlated systems. 
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